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We present a quantization scheme for the electromagnetic field interacting with atomic systems 
in the presence of dispersing and absorbing magnetodielectric media, including left-handed material 
having negative real part of the refractive index. The theory is applied to the spontaneous decay 
of a two-level atom at the center of a spherical free-space cavity surrounded by magnetodielectric 
matter of overlapping band-gap zones. Results for both big and small cavities are presented, and 
the problem of local-field corrections within the real-cavity model is addressed. 
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I. INTRODUCTION 



The problem of propagation of electromagnetic waves 
in materials having, in a certain frequency range, simul- 
taneously negative permittivity and permeability thus 
leading to a negative refractive index was first studied 
by Veselago Q . Since in such materials the electric field, 
the magnetic field, and wave vector of a plane wave form 
a left-handed system, so that the direction of the Poynt- 
ing vector and the wave vector have opposite directions, 
they are also called left-handed materials (LHMs). Other 
unusual properties are a reverse Doppler shift, reverse 
Cerenkov radiation, negative refraction, and reverse light 
pressure. Since LHMs do not exist naturally, they have 
remained a merely academic curiosity until recent reports 
on their fabrication [3- C3- 123- 13- IZ3| - The metamaterials con- 
sidered there consist of periodic arrays of metallic thin 
wires to attain negative permittivity, interspersed with 
split-ring resonators to attain negative permeability. Al- 
though the metamaterials that have been available so 
far behave like LHMs only in the microwave range, there 
have been suggestions on how to construct metamaterials 
that can operate at optical frequencies, by reducing the 
sizes of the inclusions (split rings, chiral or omega par- 
ticles) 7] or by using point defects in photonic crystals 
as magnetic emitters ||. A number of potential appli- 
cations of LHMs have been proposed, including effective 
light-emitting devices, beam guiders, filters, and near- 
field lenses. For example, LHMs could be used to real- 
ize highly efficient low reflectance surfaces |jj or super- 
lenses which, in principle, can achieve arbitrary subwave- 
length resolution |l(j . The intriguing superlense proposal 
and the reported observation of negative refraction |3J 
have touched off intensive and enlightening discussions 



[Til IT3 . Illij ] . More recent experiments @ seem to confirm 
the negative refraction observed in Ref. 3]. Neverthe- 
less, there have been still many open questions about 
the electrodynamics in magnetodielectrics, i.e., materi- 
als with simultaneously significant electric and magnetic 
properties, including LHMs. 

In this paper, we first study the problem of quanti- 
zation of the macroscopic electromagnetic field in the 
presence of magnetodielectrics, with special emphasis on 
LHMs. Apart from the more fundamental interest in the 
problem, quantization is required to include nonclassical 
radiation in the studies. Since dispersion and absorp- 
tion are related to each other by the Kramers-Kronig 
relations, noticeable dispersion implies that absorption 
also cannot be omitted in general. As we will show, 
quantization of the electromagnetic field in the presence 
of dispersing and absorbing magnetodielectrics can be 
performed by means of a source-quantity representation 
based on the classical Green tensor in a similar way as in 
Refs. [HE [HE El for P urel y dielectric material. 

As a simple application of the quantization scheme, we 
then study the spontaneous decay of an excited two-level 
atom in a dispersing and absorbing magnetodielectric en- 
vironment, with special emphasis on an atom in a spher- 
ical cavity. It is well-known that the spontaneous decay 
of an atom is influenced by the environment. If the atom 
is embedded in a homogeneous, purely electric medium 
with real and positive (frequency-independent) permit- 
tivity, the decay rate without local-field corrections reads 



L = nLn 



(1) 
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where Fo is the decay rate in free space and n = y/e is the 
refractive index (see, e.g., f20t f2l| and references therein). 
From energy scaling arguments it can be inferred that the 
electric field in a medium corresponds to the electric field 
in free space times 1 / ^fe. From a mode decomposition 
one can conclude that the mode density is proportional to 
n 3 . With that, Eq. immediately follows from Fermi's 
golden rule. Now if we take into account that in the 



2 



more general case of positive e and [i the refractive index 
is n = y/sjl, we conclude that 



r = finTo- 



(2) 



Unfortunately, these arguments cannot be used if, e.g., 
fi and n are simultaneously negative. Basing the calcu- 
lations on rigorous quantization, we show that Eq. PJl 
also remains valid in this case. Moreover, we generalize 
Eq. (J2J) to the realistic case of dispersing and absorbing 
matter, including local-field effects. 

The article is organized as follows. In Sec.[nJ some gen- 
eral aspects of the refractive index of a medium whose 
permittivity and permeability can simultaneously be- 
come negative are discussed. Section ITTT1 is devoted to the 
quantization of the electromagnetic field in the presence 
of a dispersing and absorbing magnetodielectric medium. 
The interaction of the medium-assisted field with addi- 
tional charged particles is considered in Sec. HVI and the 
minimal-coupling Hamiltonian is given. In Sec. the 
theory is applied to the spontaneous decay of an excited 
two-level atom, with special emphasis on an atom in a 
spherical cavity surrounded by a dispersing and absorb- 
ing magnetodielectric. A summary and some concluding 
remarks are given in Sec. IV1I 



II. PERMITTIVITY, PERMEABILITY, AND 
REFRACTIVE INDEX 



Let us consider a causal linear magnetodielectric me- 
dium characterized by a (relative) permittivity er(r, uj) 
and a (relative) permeability fi(r, uj), both of which are 
spatially varying, complex functions of frequency satisfy- 
ing the relations 

e(r, -uj*) = e*(r,cj), /j(r, -u*) = //*(r, uj). (3) 

They are holomorphic in the upper complex half-plane 
without zeros and approach unity as the frequency goes 
to infinity, 

lim e(r, u>) = lim fi(r, uj) = 1. (4) 

\u\ — >oo \u\ — >oo 

Since for absorbing media lme(r, uj) > 0, Im/i(r, uj) > 
(see, e.g., Ref. H3), we may write 

e(r,uj) = \e(r,uj)\e^^\ <f> e (r, uj) G (0, tt), (5) 
f,(r,uj) = \fi(r,uj)\e i ^ u \ ^(r, uj) G (0, tt). (6) 

The relation n 2 (r,uj) =e(r, uj)fi(r, uj) formally offers two 
possibilities for the (complex) refractive index n(r, uj), 



n(r,w) = ±y/\e(r, wjJ^F, cu)\ jlM'rf+MwW, (7) 
where 

O<[0 e (r,«)+^(r,o;)]/2<7r. (8) 



The ± sign in Eq. (J7J) leads to Imu(r,w) ^ 0. To specify 
the sign, different arguments can be used. From the high- 
frequency limit of the permittivity and the permeability, 
Eq. @, and the requirement that limi^i—joo n(r, uj) = 1 it 
follows that the + sign is correct |l2(, 

n(r,uj) = y/\e(r, ufifir, w)\ e »[<M r ^)+"M^)]/2, (9) 

The same result can be found from energy arguments [t| 
[see also the remark following Eq. i|27|) in Sec. II I II and 
Sec. ESI- 

From Eq. @ it can be seen that when both e(r,uj) 
and /i(r, uj) have negative real parts [4> e {r,uj), ^(r, uj) G 
(7r/2,7r)], then Ren(r, w) is also negative. It should be 
pointed out that for negative Re n(r, uj) it is not nec- 
essary that Ree(r, uj) and Re^i(r,w) are simultaneously 
negative. For the real part of the refractive index to 
be negative, it is sufficient that [<f> e (r,uj) + p (r, uj)} > 
tt, i.e., one of the phases can still be smaller than 7r/2, 
provided the other one is large enough. In fact, the def- 
inition of LHMs was originally introduced for frequency 
ranges where material absorption is negligible small, and 
thus e(r,uj) and n{r,uj) can be regarded as being real 
In this case propagating waves can exist provided 
that both e(r, uj) and /i(r,uj) are simultaneously either 
positive or negative. If they have different signs, then 
the refractive index is purely imaginary, and only evanes- 
cent waves are supported. The situation becomes more 
complicated when material absorption cannot be disre- 
garded, because there is always a nonvanishing real part 
of the refractive index (apart from the specific case where 
[4> e (r,uj) + <j>fj,(r,uj)] =7r). In the following we refer to a 
material as being left-handed if the real part of its refrac- 
tive index is negative. 

In order to illustrate the dependence on frequency of 
the refractive index, let us restrict our attention to a 
single-resonance permittivity 



e{u) = I + 



and a single-resonance permeability 



i 



(10) 



(11) 



where ujp e , uip m are the coupling strengths, uj^ e , uj^m 
are the transverse resonance frequencies, and j e , j m are 
the absorption parameters. For notational convenience, 
we have omitted the spatial argument. Both the permit- 
tivity and the permeability satisfy the Kramers-Kronig 
relations. Equation (|10|l corresponds to the well-known 
(single-resonance) Drude-Lorentz model of the permittiv- 
ity. The permeability given by Eq. lllfl . which is of the 
same type as Eq. (|10f) . can be derived by using a damped- 
harmonic oscillator model for the magnetization [23|. It 
also occurs in the magnetic metamaterials constructed 
recently HSHEI. 

For very small material absorption (j e / m <C uj P e / m , 
w Te/m)j the permittivity l|l(J|l and the permeability 
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FIG. 1: Real (a) and imaginary (b) parts of the refractive 
index n(ui) as functions of frequency, with the permittivity 
e{uj) and the permeability /j.(uj) being respectively given by 

EqS. itTOl and {HJ [u>Te = l.03iUTm; Wp m = 0.43 LdTm', UPe = 

0.75 WTm! — — 0.001 uJTm (solid lines) , 0.01 uJTm (dashed 
lines), and 0.05 u>Tm (dotted lines)]. The values of the param- 
eters have been chosen to be similar to those in Refs. la. Ha. 



respectively, feature band gaps between the transverse 
frequency WTe and the longitudinal frequency WL e = 



\J uj\ e + w Pe' wnere Ree(w) <0, and the transverse fre- 
quency LOTm and the longitudinal frequency WL m = 
V^Tm + W L> where Re h(uj) < 0. With increasing val- 
ues of 7e and 7 m the band gaps are shifted to higher 
frequencies and smoothed out. Figure [I] shows the de- 
pendence on frequency of the refractive index [Eq. J3J] 
for the case of overlapping band gaps and various ab- 
sorption parameters. In particular, if max (a>i>e, w Tm) < 
lo < min (c^Le, WLm), then Ree^) < and Re/i(cj)<0, 
and thus a negative real part of the refractive index is 
observed. In Fig. this is the case in the frequency 
interval where 1.03 < w/wTm < 1.088. For the chosen 
parameters, a negative real part of the refractive index 
can also be realized for frequencies slightly smaller than 
ujT e , where Refi(uj) < while Ree(uj) > 0, as is clearly 
seen from the inset in Fig. Qla). In this region, how- 
ever, |Ren(o;)| is typically small whereas Imn(w) is large 
thereby effectively inhibiting traveling waves. It is worth 



noting that a negative real part of the refractive index 
is typically observed together with strong dispersion, so 
that absorption cannot be disregarded in general. On the 
other hand, increasing absorption smooths the frequency 
response of the refractive index thereby making negative 
values of the refractive index less pronounced. 



III. THE QUANTIZED MEDIUM-ASSISTED 
ELECTROMAGNETIC FIELD 

The quantization of the electromagnetic field in a 
causal linear magnetodielectric medium characterized by 
both e(r,u>) and /j,(t,uj) can be performed by generaliz- 
ing the theory given in Refs. [la.ll9| for dielectric media. 
Let P(r, uj) and M(r,Lj), respectively, be the operators 
of the polarization and the magnetization in frequency 
space. The operator-valued Maxwell equations in fre- 
quency space then read 

VB(r,u/)=0, (12) 

VD(r,w) = 0, (13) 

V x E(r, to) = wB(r, uj), (14) 

V x H(r, uj) = -iujD(r, uj), (15) 

where 

D(r,w) = £ E(r,w) +P(r,w), (16) 

H(r,w) = KoB(r,w)-M(r,w) (17) 
[kq = Hq 1 ]- Similarly to the electric constitutive relation, 

P(r, uj) = e [e(r, uj) - l]E(r, uj) + P N (r, lu), (18) 

with P N (r,cj) being the noise polarization associated 
with the electric losses due to material absorption, we 
introduce the magnetic constitutive relation 

M(r, uj) =ko[1- w)]B(r, uj) + M N (r, uj), (19) 

where ft(r,a-i) = ^ _1 (r,u;), and M N (r, uj) is the noise mag- 
netization unavoidably associated with magnetic losses. 
Recall that for absorbing media Im/i(r, uj) > 0, and thus 
Imft(r,w) <0. Substituting Eqs. {n} and JT^-JT^ into 
Eq. 115|) . we obtain 



Vx K (r,o;)VxE(r,w) ^(r, w)E(r, uj) = iw/x j N (r, w), 

(20) 

where 



j*r( r > w ) = -iwPjvj^w) + V xM N (r,w) 



(21) 



is the noise current. The noise charge density is given 
by p (r, uj) = — VP N (r, w), and the continuity equation 
holds. The solution of Eq. (|20J) can be given by 

E(r,u;)=^Mo f dV G(r, r', to) L(r', uj), (22) 
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where G(r, r', uj) is the (classical) Green tensor satisfying 
the equation 



V x k(t, uj) V x e(r, ui) 



G(r,r',u))=S(r-r') 



(23) 

together with the boundary condition at infinity. It 
is not difficult to prove that the relation G*(r, r',u>) 
= G(r, r', —uj*), which is analogous to the relations J3J), 
is valid. Other useful relations are (see Appendix 



Gy(r,r» = Gji(r',r,u) 



(24) 



and 



d 3 s|-IniK(s, uj) G(r,s,w)xV s V s x G*(s, r', w) 
^Im£(s,w)G(r,s,w)G , (s 1 r»} = ImG(r,r' 



where 



(25) 



G(r, s,w)xV s = e jk id s k G u (r, s, w). (26) 



In the simplest case of bulk material, Eq. (|23|l implies 
that the Green tensor can simply be obtaine d by multi- 
plying the Green tensor for a bulk dielectric [13, Wl E3 
by /x(w) and replacing e(w) with s(uj)fj,(u>), 



G^(r,r» =h(lj) [flTflf + g 2 M%(r - r')] 



giRe g(a;)|r— r' 1 

47rg 2 (^)|r - r'| 



-Im q(io) |r— r \ 



(27) 



[q(uj) = n{w)u>/c\. From the boundary condition for the 
Green tensor at |r — r'| — > 00, it follows that Imn(u) > 0, 
which is consistent with Eq. |JS}. 

Analogously to the noise polarization that can be re- 
lated to a bosonic vector field f e (r, uj) via 



P N (r,w)=ty^Ime(r,«)t s (r > a;), (28) 

the noise magnetization can be related to a bosonic vector 
field f m (r, w) via 



M N (r,w)= l /-^ImKir„.)f,„( r ...]. (2!)) 
with (A, A' = e, m) 

/ Ai (r,cj)jt,.(r>')l =5 A vM(r-r')<5(w-o;'), (30) 



f Xi (r,uj),fx'j(r',uj') 



= 0. 



(31) 



Substituting in Eq. ij?T]l for P N (r,w) and M N (r, uj) the 
expressions (I28f) and (|29|) . respectively, we may express 



j N (r, uj) in terms of the bosonic fields f\(r, uj) as follows: 



j x ( r. i = ^ Ime(r, w) f e (r, w) 



+ V x \J-^lmK(r,uj)i m (r,uj). (32) 

Note that in Eqs. and respectively, P N (r,w) 

and M N (r, w) are only determined up to some phase fac- 
tors which can be chosen independently of each other. 
Here we have them chosen such that in Eq. 132|) the co- 
efficients of f e (r,w) and f m (r, uj) are real. 

The f\(r,uj) and fj(r, uj) can be regarded as being 
the fundamental variables of the system composed of the 
electromagnetic field and the medium including the dis- 
sipative system, so that the Hamiltonian can be given by 



H= J ' ^ r J dwfiwfJ(r,w)f A (r,w) 



(33) 



In this approach, the medium-assisted electromagnetic- 
field is fully expressed in terms of the f A (r, uj) and 
f\(r, uj). In particular, the electric-field operator (in the 
Schrodinger picture) reads 



E(r) 



dwE(r,w) + H.c. 



(34) 



where E(r, uj) is given by Eq. (|2*2"|) together with Eq. 
Similarly, the other fields can be expressed in terms of the 
f\(r, uj) and fJ(r,o;), by making use of Eqs. I|14|l. IjlSjl. 
$P$ , iPBJ, and It can then be shown (Appendix 

IB1 that the fundamental (equal-time) commutation rela- 
tions 

[E i (r) i E j (r')]=Q=[B i (r),B j (r')], (35) 
[e Ei(r), 4-(r')] = -ihe ijk d r k 6(v - r') (36) 

are preserved. Furthermore, it can be verified (Appendix 
[UJ) that in the Heisenberg picture the medium-assisted 
electromagnetic-field operators obey the correct time- 
dependent Maxwell equations. 

The introduction of a noise magnetization of the type 
of Eq. H29|) was first suggested in Ref. but it was 
wrongly concluded that such a noise magnetization and 
the noise polarization in Eq. (|28|l can be related to a 
common bosonic vector field f(r, uj). Since f e (r,o;) in 
Eq. (|28(l is an ordinary vector field, whereas f m (r,cj) in 
Eq. 129|) is a pseudo-vector field, the use of a common 
vector field would require a relation for the noise magne- 
tization that is different from Eq. (|29|l but must ensure 
preservation of the commutation relations i|35[) and i|36|) 
and lead to the correct Heisenberg equations of motion. 
For the metamaterial considered in Refs. 0, 0, 0, IE 0- 
where the electric properties and the magnetic properties 
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are provided by physically different material components, 
the assumption that the polarization and the magnetiza- 
tion are related to different basic variables is justified. It 
is also in the spirit of Ref. |2^|, where the polarization 
and the magnetization are caused by different degrees of 
freedom. 



IV. INTERACTION OF THE 
MEDIUM-ASSISTED FIELD WITH CHARGED 
PARTICLES 

In order to study the interaction of charged particles 
with the medium-assisted electromagnetic field, we first 
introduce the scalar potential 



the third term is their Coulomb energy, with 



0(r) = J du> <p(r, lj) + H.c. 



and the vector potential 



A(r) = / dwA(r,w) +H.c, 



(37) 



(38) 



where in the Coulomb gauge, </3(r, ui) and A(r, u)) are 
respectively related to the longitudinal part E' (r, uj) and 
the transverse part E (r, w) ofE(r, u>) [Eq. ll'L'l) together 
with Eq. glj] according to 



-Vc/3(r,w) =E ll (r,w), 
A(r,cj) = (iw)" 1 E ± (r,w). 



(39) 
(40) 



Similarly, the momentum field II(r) that is canonically 
conjugated with respect to the vector potential A(r) can 
be constructed noting that n(r, u>) = — £olf L ( r > w). Now 
the Hamiltonian (|33|l can be supplemented by terms de- 
scribing the energy of the charged particles and their in- 
teraction energy with the medium-assisted electromag- 
netic field in the same way as in Ref. for dielectric 
matter. In the minimal-coupling scheme and for non- 
relativistic particles, the total Hamiltonian then reads 

H = / d3r / duJ HuJ *A ( r > u ^ ( r ' u ) 



A— e,m 



E 



2 m n 



+ iy d 3 rp A (r)£ A (r) + J d 3 rp A (r)<£(r), (41) 

where f Q , and p a are respectively the position and the 
canonical momentum operator of the ath particle of mass 
m a and charge q a . The first term in Eq. I|41(l is the 
Hamiltonian (|33J) of the electromagnetic field and the 
medium including the dissipative system. The second 
term is the kinetic energy of the charged particles, and 



Mr) = ^2q a 6(r-r a ), 



d 3 r / PA(r') 



47T£n r — r' 



(42) 
(43) 



being respectively the charge density and the scalar po- 
tential of the particles. Finally, the last term is the 
Coulomb energy of the interaction between the charged 
particles and the medium. 

Let E(y) and B(r) be respectively the operators of the 
electric field and the induction field in the presence of the 
charged particles 



E(r) - E(r) - V^ A (r), B(r) = B(r) 



(44) 



Accordingly, the displacement field D(r) and the mag- 
netic field H(r) in the presence of the charged particles 
are given by 

B(t) = D(r) - e V^ A (r), l(r) - H(r). (45) 

Note that in Eqs. (|44|> and (|45|l the electromagnetic fields 
must be thought of as being expressed in terms of the fun- 
damental fields f\(r) and ft(r). From the construction 
of the induction field and the displacement field it follows 
that they obey the time-independent Maxwell equations 



VB(r) = 0, VD(r) = p A (r) 



(46) 



Further, it can be shown (Appendix Q that the Hamil- 
tonian 1)41(1 generates the correct Heisenberg equations of 
motion, i.e., the time-dependent Maxwell equations 



V x E(t)+B(t) = 0, 
Vxl(r)J(r)=j A (r), 



(47) 
(48) 



where 



j A (r) = \ q a [r a 5(r -r a ) + S(r- r a )h a , (49) 



and the Newtonian equations of motion for the charged 
particles 



1 



rn 



IV - q a A(i a ) 



(50) 



a = q a (E(r a ) + \ r a x B(r a ) - B(r a ) x f „ | 



(51) 
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V. SPONTANEOUS DECAY OF AN EXCITED 
TWO-LEVEL ATOM 



(|56ll and substitution into Eq. i|54|) leads to, upon using 
the relation Q25p. 



Let us consider a two- level atom (position r^, tran- 
sition frequency u>a) that resonantly interacts with 
the electromagnetic field in the presence of magnetodi- 
electrics and restrict our attention to the electric-dipole 
and the rotating-wave approximations. By analogy with 
the case of an atom in the presence of dielectric material 
[T^,l26l| . the Hamiltonian (|41f) reduces to 



H 



E 

A— e,m 



d 6 r 



<7 f cL 



duj Huj t' x (r, w)f\ (r, uj) 



duE(r A ,w) + H.c. 



(52) 



where a = \l){u\ and cr = \u)(l\ are the Pauli operators 
of the two-level atom. Here, |Z) is the lower state whose 
energy is set equal to zero and \u) is the upper state 
of energy fiwA- Further, <1a = (/|c!a|u) = (u|c1a|0 is the 
transition dipole moment. 

To study the spontaneous decay of an initially excited 
atom, we may look for the system wave function at time 
t in the form of [|l A (r, u)) = fj>, u/)|{0})] 

m))=Cu(t)e-^ t \{0})\u) 

+ E / d3r / dwe-^CAjfow.tJIlAfow))!*), (53) 

A— e,m 

where C u (t) and C\i(t) are slowly varying amplitudes 
and, in anticipation of the environment-induced tran- 
sition frequency shift Sui |27|. uja = uja — Suj is the 
shifted transition frequency. The Schrodinger equation 
ifodt\4>(t)} = H\ip(t)) then leads to the set of differential 
equations 



C u {t) = -iSujC u (t) 



J d 3 r d A | ^ ^lm.e(r, w) G(r A , r, uj)C e i(r, u, t) 
y/-ImK,(r,u>) G(r A ,r,w)xV r C ml (r,u,t) >, (54) 



C e l(r,LU,t) 



C m i(r,w,t) 



1 uj 2 



xd A G*(v A ,r,uj)C u (t), 
y/—Im k(t, uj) e 



(55) 



\/Tre h c 



i(u—a>A)t 



x d/s 



G*(r A ,r,^) x V r C u (t), (56) 



which has to be solved under the initial conditions C u (0) 
= 1, CA;(r,w,0) = 0. Formal integrations of Eqs. I|55|l and 



C u (t) = -i5ujC u {t) + / dt' K(t - t')C u {t'), (57) 
Jo 



where 



K{t-t') 



u> 



h-KEo J C 2 



x e-^—^JCt-f ) dAlm G(rAj rAj w)dA , ( 58) 

It should be noted that, by integrating with respect to t, 
the integro-differential equation (|57|l can equivalently be 
expressed in the form of a Volterra integral equation of 
second kind [2(|. Equations (|57|l and (|58|l formally look 
like those valid for non-magnetic structures [2(j. Since 
the matter properties are fully included in the Green ten- 
sor, the results only differ in the actual Green tensor. 

Equations 157|l and (|58() apply to an arbitrary cou- 
pling regime [26J. Here, we restrict our attention to the 
weak-coupling regime, where the Markov approximation 
applies. That is to say, we may replace C u (t') in Eq. (|57|) 
by C u (t) and approximate the time integral according to 



Jo 



(59) 



[C(x) =tt5(x) +iP/x\. Identifying the principal-part in- 
tegral with the transition- frequency shift, we obtain 



Slu = 



1 



nhe. 



o Jo 



u) 2 d A IniG(rA,rA,w)d/ 



UJ — U>A 



, (60) 



which, together with Eq. I|53|). can be regarded as being 
the self-consistent defining equation for the transition- 
frequency shift 27] . Equation (|57|l then yields C u (t) — 
exp(— |rt), where the decay rate T is given by the formula 



2ujf_ 
Heoc 2 



d A Im G(r A , r A , w A )d A , 



(61) 



which is obviously valid independently of the (material) 
surroundings of the atom. 



A. Nonabsorbing bulk material 

Let us first consider the limiting case of non-absorbing 
bulk material, i.e., s(uj a ) and h(uj a ) are assumed to be 
real. Using the bulk-material Green tensor (|27|l . it can 
easily be proved that 

Im G(r A , r A , uj a ) = Re [n(uj A )n(u) A )] I. (62) 

67TC 

Substitution into Eq. l|6*T|) yields the decay rate 

r = Re[/i(£ A Mw A )]r , (63) 
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where 



r = 



3hireoc 3 



(64) 



is the free-space decay rate, but taken at the shifted tran- 
sition frequency. Eq. I|63|) is in agreement with Eq. (J2J) ob- 
tained from simple arguments on the change of the energy 
density and the mode density for positive and frequency- 
independent e and /i. Clearly, Eq. (I<j3|> is more general in 
that it also applies to dispersive magnetodielectrics. In 
particular, when e(Qa) and //((Da) have opposite signs, 
then the refractive index defined according to Eq. © 
is purely imaginary, thereby leading to T = 0. This is 
because the electromagnetic field cannot be excited at 
d>A, so that spontaneous emission is completely inhib- 
ited. Note that material absorption always gives rise to 
a finite value of T, which of course can be very small. 

From Eq. (|63[) it is clearly seen that for non-absorbing 
LHM, i.e., e(<Da) < and /x(u)a) < 0, the now real refrac- 
tive index must also be negative, in order to arrive at a 
non-negative value of the decay rate. This is yet another 
strong argument for the choice of the + sign in Eq. ■ 



B. Atom in a spherical cavity 



For realistic bulk material, the imaginary part of the 
Green tensor at equal positions is singular 0, 0, |2{| . 
Physically, this singularity is fictitious, because the atom, 
though surrounded by matter, is always localized in a 
small free-space region. The Green tensor for such an 
inhomogeneous system reads 

G(r, r',u>) = G v (r, r', u) + G s (r, r', u>), (65) 

where G v (r, r\ui) is the vacuum Green tensor and 
G s (r,r',w) is the scattering part, which describes the 
effect of reflections at the surface of discontinuity. Using 
Eq. jnni together with ImG v (r, r,w) = (u/6irc)I [cf. 
Eq. 162|) ]. we can write the decay rate (|61|l as 



r = r n 



heoc 2 



d A Im G s (r A , r A , <Da)cIa, 



(66) 



which is again seen to be valid for any type of material. 

Within a 'classical' theory of spontaneous emission 
[28| , a formula of the type (|66[) has been used in Ref . [2j| 



to calculate the decay rate of an atom near a dispersion- 
less and absorptionless LHM sphere. 'Classical' theory 
means here, that a classically moving dipole in the pres- 
ence of macroscopic bodies is considered, with the value 
of To being borrowed from quantum mechanics. As in 
Ref. pflj . the atomic transition frequency is commonly 
understood as being that in free space. From Eq. 1(66(1 it 
is seen that the medium-assisted (i.e., shifted) frequency 
(Da must be used instead of the free-space frequency wa, 
since both values can differ substantially. 

Let us apply Eq. 1)66(1 to an atom in a free-space re- 
gion surrounded by a multilayer sphere. Using the Green 
tensor given in Ref. [30| , we obtain 



= l + f^n(rc+l)(2n+l) 



-i 2 



^TA 



ReC 
(67) 



for a radially oriented dipole moment (dA || i"a) and 



- = l + |£(2n + l) 



^(/c A r A )ReC, 



M 



[k A rAjn(k A rA)}' 

k A TA 



Re C, 



N 



(68) 



for a tangentially oriented dipole moment (dA -Lta) [the 
prime indicating the derivative with respect to k A r A , 
(fc A = w A /c)]. In Eqs. ED and JBSJ, j n (z) and h£\z) 
are the spherical Bessel and Hankel functions of the first 
kind, respectively. The coefficients and Cjf have to 
be determined through recurrence formulas |30j . 

Equations (|67|l and l|68|) apply to an atom at an ar- 
bitrary position inside a spherical free-space cavity sur- 
rounded by an arbitrary spherical multilayer material en- 
vironment. Let us specify the system such that the atom 
is situated at the center of the cavity (i.e., ta = 0) and 
let the surrounding material homogeneously extend over 
all the remaining space. For small cavity radii, the sys- 
tem corresponds to the real-cavity model of local-field 
corrections. Making use of the explicit expressions for 
the coefficients as in Ref. [3(j and the fact that for 
ta = only the n = 1 term in Eq. I|67|) contributes [3l| , 
we derive from Eq. (|67|l 



T 



1 + Re 



1 - i(n+ l)z - n(n + 1) ^ " ' z 2 + in 2 - - - 3 



fi—n 



-isinz — (nsinz — icosz)z + ( cos z — i ^-nsinz ] nz 2 — (nsinz + iucosz) -z 3 

fi — n 2 J [i — n 2 - 



(69) 



\jjl = /i(cja), n = n{QA)i and z = Ruja/c, with R being the cavity radius]. Obviously, the dipole orientation does not 
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FIG. 2: The decay rate T as a function of the (shifted) atomic 
transition frequency o>a for an atom at the center of an empty 
sphere surrounded by single-resonance matter, (a) dielectric 
matter according to Eq. (| 1 ( ) p [wTe/wTm = 1.03; wpe/tWm = 
0.75; 7 e /ti)Tm = 0.001 (solid line), 0.01 (dashed line), and 
0.05 (dotted line)], (b) magnetic matter according to Eq. pip 
[ujp m /ujTm = 0.43; 7 m /u>Tm = 0.001 (solid line), 0.01 (dashed 
line), and 0.05 (dotted line)], and (c) magnetodielectric mat- 
ter according to Eqs. 11 0ft and pip [the parameters are the 
same as in (a) and (b)] . The diameter of the sphere is 2R — 

20 Aim (XTm — 2-Kc/iJTm)- 



matter here, and Eqs. (|67j) and I|fi8|l lead to exactly the 
same result. Equation 169(1 is the generalization of the 
result derived in Ref. |3l| for dielectric matter. 

Figures El - El illustrate the dependence of the decay 
rate T given by Eq. l|o9[) on the (shifted) transition fre- 
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FIG. 3: The decay rate T as a function of the (shifted) atomic 
transition frequency o>a for an atom at the center of an empty 
sphere surrounded by single-resonance matter, (a) dielectric 
matter according to Eq. il l (b) magnetic matter according 
to Eq. itm . and (c) magnetodielectric matter according to 
Eqs. pop and pip ['y e /^Tm = 7m/<^Tm = 0.001, the other 
parameters are the same as in Fig. H]. The values of the 
sphere diameter are 2R = 20 Ax m (dashed lines) and 2R = 
1 At™ (solid lines). 



quency for the case of the cavity being surrounded by 
(a) purely dielectric matter, (b) purely magnetic mat- 
ter, and (c) magnetodielectric matter near the band-gap 
zones. The permittivity and permeability are given by 
Eqs. (fl"n)l and fTT|l . respectively. In the figures, the di- 
electric and magnetic band gaps are assumed to extend 
from <x>Te = 1-03 o>Tm to WLe — l-274o;Tm and from LUTm 
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to WLm — 1-088 WTm, respectively. They overlap in the 
frequency interval 1.03 wxm < w < 1-088 wxm- 



1. Large cavities 

In Fig. |2 a relatively large cavity is considered 
(2R/X Tm = 20). From Figs.^a) and^b) it is seen that 
inside a dielectric or magnetic band gap the decay rate 
sensitively depends on the transition frequency. Narrow- 
band enhancement of the spontaneous decay (r > Tq) al- 
ternates with broadband inhibition (r<ro). The max- 
ima of enhancement are observed at the frequencies of 
the (propagating-wave) cavity resonances, the Q factors 
of which are essentially determined by the material losses 
(see the curves for different values of 7 e and j m ). Note 
that the cavity resonances as the poles of T are different 
for dielectric and magnetic material in general. From 
Fig- GI C ) it is seen that the decay rate of an atom sur- 
rounded by magnetodielectric matter shows a similar be- 
haviour as in Figs.|2ta) and [2(b), provided that the tran- 
sition frequency is outside the region of overlapping di- 
electric and magnetic band-gap zones. When the transi- 
tion frequency is in the overlapping region of the gaps, 
then the medium becomes left-handed. Thus, a relatively 
large input-output coupling due to propagating waves in 
the medium becomes possible, thereby the typical band- 
gap properties getting lost. As a result, neither strong 
inhibition, nor substantial resonant enhancement of the 
spontaneous decay is observed, as is clearly seen from 
Fig. He). 

In Fig. |3 the results for the cavity in Fig. [2 are com- 
pared with those observed for a smaller cavity with 
2R/X^ m = 1. As expected, the number of clear-cut cavity 
resonances decreases as the radius of the cavity decreases. 
For the smaller of the chosen radii, just one resonance has 
survived in the case of the magnetic medium [Fig. [3(b)], 
while the resonances are gone altogether in the case of the 
dielectric medium [Fig. [3(a)]. Accordingly, inhibition of 
spontaneous decay is typically observed in the band-gap 
zones of dielectric and magnetic matter and in the non- 
overlapping band-gap region of magnetodielectric matter. 
In contrast, a behaviour quite similar to that in free space 
can be observed in the overlapping (left-handed) region. 



2. Small cavities 
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FIG. 4: The decay rate T as a function of the (shifted) atomic 
transition frequency u>a for an atom at the center of an empty 
sphere surrounded by single-resonance matter, (a) dielectric 
matter according to Eq. HUH , (b) magnetic matter accord- 
ing to Eq. 1111 . and (c) magnetodielectric matter according 
to Eqs. IIUI and 1111 . The diameter of the sphere is 2R = 
0.1 At™- The other parameters are the same as in Fig. [5] 



In Fig. 0J a cavity is considered whose radius is much 
smaller than the transition wavelength (2i?/Axm = 0.1). 
Comparing Fig. 0(a) with0(b), we see that the frequency 
response of the decay rate in the dielectric band-gap zone 
is quite different from that in the magnetic band-gap 
zone. In the dielectric band-gap zone [Fig. Ufa)], a more 
or less abrupt decrease of T below Tq with increasing 
transition frequency is followed by an increase of T to a 
maximum that can substantially exceed Tq. In the case 
of magnetic matter, [Fig. 0(b)], on the contrary, only a 



rather distorted band-gap zone is observed in which T 
monotonously decreases below Tq. The maximum of en- 
hancement of spontaneous decay in Fig. 0(a) is observed 
at the local-mode resonance associated with the small 
cavity, which may be regarded as being a defect of the 
otherwise homogeneous dielectric. This is obviously of 
the same nature as the donor and acceptor local modes 
discussed in Ref. ■ In the regions where the dielectric 
and magnetic band-gap zones of the magnetodielectric in 
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FIG. 5: The decay rate V as a function of the (shifted) 
atomic transition frequency wa for an atom at the center of 
an empty sphere surrounded by single-resonance matter, (a) 
dielectric matter according to Eq. (IIOP . (b) magnetic matter 
according to Eq. 1111 . and (c) magnetodielectric matter ac- 
cording to Eqs. (1101 and Hill [j e /ujTm = 7m/wTm = 0.001, 
the other parameters are the same as in Fig.lU. The values of 
the sphere diameter are 71 = 0.8 At™ (dotted lines), 0.4 At™ 
(dashed lines), and 0.1 At™ (solid lines). 



Fig. 0fc) do not overlap, the frequency response of the 
decay rate is dominated by the respective matter, i.e., the 
characteristic features are either dielectric or magnetic. 
The situation changes when the transition frequency is 
in the overlapping region, where LHM is realized. Since 
this region cannot longer be regarded as an effectively 
forbidden zone for propagating waves, the value of T can 



become comparable with or even bigger than that of Tq. 
From Fig. 0Jc) it is seen that entering the overlapping 
region from the magnetic side stops the decrease of T on 
that side, thereby changing it to an increase. Similarly, 
the decrease of T on the dielectric side stops and changes 
to an increase when the overlapping region is entered 
from the dielectric side. 

Figure[S]iriustrates the influence of the cavity radius on 
the decay rate for small cavities. FigureEt a ) reveals that 
when the value of 2i?/Ax m changes from 2i?/ATm = 0.1 
to 2Rj Axm = 0.8, then the maximum of the spontaneous 
decay rate associated with the local-mode resonance in 
dielectric matter shifts towards smaller transition fre- 
quencies, thereby being reduced. In the case of magnetic 
matter, increasing value of 2R/\T m reduces the distor- 
tion of the band-gap zone, as is seen from Fig. Et^) . As 
expected, the frequency response of the decay rate shown 
in Fig. Etc) for the case of magnetodielectric material in- 
cluding LHM combines, in a sense, the respective curves 
in Figs.Et a ) andEtb). 



3. Local-field corrections 

For an atom in bulk material, the local field with which 
the atom really interacts can differ from the macroscopic 
field used in the derivation of the decay rate of the form 
given by Eq. (|63|) . To include local-field corrections in 
the rate, one can use Eq. (|69[) and let the radius of the 
cavity tend to a value which is much smaller than the 
transition wavelength, 



Rui A _ 2nR 
c A A 



« 1, 



(70) 



but still much larger than the distances between the 
medium constituents to ensure that the macroscopic the- 
ory applies. In this way we arrive at the real-cavity model 
frequently used in the literature S EB] The 
results shown in Fig. 0] may be regarded as being typical 
of the real-cavity model. 

Expanding T [Eq. i|fci9|) ] in powers of z = RCja/c we 
obtain 



— = Re 

r 



■ Im 



3c 



l + 2e 



fin 



9Ime 



1 



e(l + 3e + 5/ue) 



(l + 2e) 2 



\CoaR 



2e\ 2 \CoaR, 

+ 0(R), (71) 



which for nonmagnetic media reduces to results obtained 
earlier |jn],[35|]. Note that the actual value of R, which 
is undetermined within the real-cavity model, should be 
taken from the experiment. Equation (|71JI without the 
O(R) term has to be employed with great care, because 
it fails when, for small absorption, the atomic transition 
frequency u>a becomes close to a medium resonance fre- 
quency such as or wxm, thus leading to a drastic 
increase of the first term in Eq. (|71|l . The first three 
terms on the right-hand side in Eq. Q71[) reproduce the 
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curves in Fig. 0] sufficiently well, except in the vicinities 
of uj^ e and CL>Tm- ^ n particular, it can easily be checked 
that the position of the local-mode-assisted maximum of 
the decay rate in the dielectric band-gap zone is where 
2e(w A )~-l. 

For transition frequencies that are sufficiently far away 
from a medium resonance frequency and, in case of di- 
electric and magnetodielectric matter, the local-mode fre- 
quency, so that material absorption can be disregarded, 
the first term in Eq. I|71l) is the leading one, hence 



F 



1 2 



1 + 2e(w A ) 



Re[MwA)n(wA)]r . (72) 



In this case, the local-field correction simply results in 
multiplying the rate obtained for the case of nonabsorb- 
ing bulk material [Eq. (§3)] by the factor [3e/(l + 2e)] 2 . 
Interestingly, this factor is exactly the same as that for 
dielectric material. 

Inspection of the second and the third term in Eq. i|71|) 
shows that such a separation is no longer possible when 
material absorption must be taken into account. It 
should be pointed out that the second term proportional 
to i?~ 3 is purely dielectric, whereas the magnetization 
starts to come into play only via the third term propor- 
tional to J? -1 . These two terms can be regarded as re- 
sulting from the near-field component and the induction- 
field component accompanying the decay of the excited 
atomic state. In particular for sufficiently small cavity 
size and strong (dielectric) absorption, the second term 
is the leading one, so that magnetodielectrics approxi- 
mately give rise to the same decay rate as dielectrics: 



r 



9Ime(d)A) 
\l + 2e(u A )\ 2 \u A R 



r . 



(73) 



In this case, the decay may be regarded as being purely 
radiationless, with the energy being transferred from the 
excited atomic state to the surrounding medium medi- 
ated by the near field. 



VI. SUMMARY AND CONCLUSIONS 

It has been shown that the quantization scheme origi- 
nally developed for the electromagnetic field in the pres- 
ence of dielectric matter described in terms of a spa- 
tiall y v arying, Kramers-Kronig-consistent permittivity 
[Tl TT5L ITa Il7l HI ITgj can be extended to causal mag- 
netodielectric matter, with special emphasis on the re- 
cently fabricated metamaterials, including LHMs that 
can exhibit a negative real part of the refractive index, 
thereby leading to a number of unusual properties. The 
quantization scheme is based on a source-quantity rep- 
resentation of the medium-assisted electromagnetic field 
in terms of the classical Green tensor and two indepen- 
dent infinite sets of appropriately chosen bosonic basis 
fields of the system that consists of the electromagnetic 



field and the medium, including a dissipative system. We 
have further shown that the minimal-coupling Hamilto- 
nian governing the interaction of the medium- assisted 
electromagnetic field with additional charged particles 
can be obtained from the standard form, by expressing 
in it the potentials in terms of the bosonic basis fields. 
The theory can serve as basis for various studies, includ- 
ing generation and propagation of nonclassical radiation 
through magnetodielectric structures, Casimir forces be- 
tween magnetodielectric bodies, or van-der-Waals force 
between atomic systems and magnetodielectric bodies. 

As an example, we have applied the theory to the prob- 
lem of the spontaneous decay of a two-level atom in the 
presence of arbitrarily configured, dispersing and absorb- 
ing media. In particular, we have shown that the theory 
naturally gives the decay rate and the frequency shift in 
terms of the classical Green tensor - formulas that are 
valid for any kind of geometry and material. To be more 
specific, we have studied the decay rate of an atom at 
the center of a cavity surrounded by a magnetodielectric, 
assuming a single-resonance permittivity and a single- 
resonance permeability of Drude-Lorentz type. LHM is 
realized for transition frequencies in the region where the 
dielectric and magnetic band-gap zones overlap, thereby 
the real parts of the permittivity and permeability be- 
coming negative. When the transition frequency enter 
that region from the dielectric or magnetic side, then the 
typical band-gap properties such as enhancement of the 
spontaneous decay at the cavity resonances and inhibi- 
tion between them gets lost and a decay rate comparable 
with that in free space can be observed. The calcula- 
tions have been performed for both large and small cavi- 
ties. In particular, if the diameter of the cavity becomes 
small compared to the transition wavelength of the atom, 
the system reduces to the real-cavity model for including 
local-field corrections in the decay rate of the atom in 
bulk material. We have discussed this case in detail both 
analytically and numerically and made contact with the 
results obtained from simple mode-decomposition argu- 
ments in case of positive permittivity and permeability. 

For simplicity, all the calculations have been performed 
for isotropic magnetodielectric material, by assuming a 
scalar permittivity and a scalar permeability. The exten- 
sion to anisotropic material is straightforward. It can be 
done in essentially the same way as for anisotropic dielec- 
tric material, by first transforming the permittivity and 
permeability tensors into their diagonal forms. 
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APPENDIX A: SOME PROPERTIES OF THE 
GREEN TENSOR 

Following Ref . 0] , we regard the Green tensor as be- 
ing the matrix elements in the position basis of a tensor- 
valued Green operator G = G(ui) in an abstract single- 
particle Hilbert space, G(r,r',ui) = (r|G|r'), so that 
Eq. (|23[l can be regarded as the position-representation 
of the operator equation HG = I, where 

uj 2 

H = p x «(f,w)p x 5-e(r,w)J. (Al) 

d- 

Using the relations (r|r|r') = r£(r — r'), (r|p|r') = 
— iV<5(r — r'), and (r| J|r') = <5(r — r'), we have 

H{r,r',uj) = (r\H\r') 

UJ 2 

= Vx K (r,w)V x 8(r-r') - _ e ( r ,wWr-r / ), (A2) 
which in Cartesian coordinates reads 
H«(r,r , ,w) = {^K(r,w)ar 



^U(r-r'). (A3) 



9 ; r K(r,w)9[ + — e(r,u>) 



Since H is injective and thus an invertible one-to-one 
map between vector functions, we can write G = H 1 . 
Multiplying this equation by H from the right, we have 



(A4) 



(A5) 



Recalling Eq. (|A3I) . we derive, on integrating by parts 
and taking into account that the Green tensor vanishes 
at infinity, 



GH = I, 

which in the position basis reads 

d 3 s(r|G|s)(s|i?|r') = 5(r - r') 



J d 3 s Gik(r, s, uj)H k j (s, r', uj) 



d[' K (r',u,)dl' + ^£(r» 



x G lk (v,v',ui) = 5ij(r - r') 



6 k j 
(A6) 



Interchanging the vector indices i and j and the spatial 
arguments r and r', we obtain 

{d r kK (r,oj)dl - [d[K(r,u)a[ + ^e{v^)]5 kl ) 

xG Jfc (r / ,r,w)=<f ii (r-r / ), (A7) 

which, according to Eq. I|23[) . is just the defining equation 
for G k j{r, r', u). Thus, the reciprocity relation l|24|) is 
proved valid. 



To prove the integral relation (|25() , we introduce oper- 
ators O* by (6*).. = (6^ = 6t. From Eq. {KQ it 



then follows that 



H X G X = I. 



(A8) 



Multiplying Eq. ljA4j) from the right by G x and Eq. (|X8|) 

from the left by G and subtracting the resulting equa- 
tions from each other, we obtain 

G(H - = G* - G, (A9) 

which in the position basis reads 

Jd 3 s J dV G im (r,s,u))[H mn (s,s',u>) 

- H* m (s',s,uj)]G* nj (s',r',u) = -2iImGy (r,r', ui). 

(A10) 

Note that <r|^ m |r') = ff* m (r',r,u;) and (r|(%|r') = 
G^r' ,r,u>). Inserting Eq. llAH into Eq. HA10|) . after 
some manipulation we derive 

d 3 s ^lmK(s,uj)d^G im (r,s,uj) 
x[d s m G* n Js,r',cj)-d^G* mj (s,r',cj)} 



uj 2 ■) 
— Im e(s, u>)G im (r, s, uj)G* mj (s, r', w) j 



ImG y (r, r',w), 



(All) 



which is just Eq. (|25|l in Cartesian coordinates. 

To examine the asymptotic behaviour of the Green ten- 
sor in the upper half of the complex w-plane as \u\ — > oo 
and \uj\ —> 0, we introduce the tensor-valued projectors 



i ± =t _P fll = E®P 



(A12) 



[note that (r|/- L (H) |r') = ^(ll)^ - r')], and decompose G 
as pL9| . 

G = H 1 = /"(/Hff/ll)- 1 /" 

+ [i± - i\\(i\\Hi\\)-ii\\Hi±]K 

x [J x - J- L fl"j"(jll J H'jll)- 1 jll] ! (A13) 



where 



K= [i^HI^ - i^H^iiWHI^^Hi 1 -} . (A14) 

Recalling that e(r,uj), /j,(r,u>) — > 1 as |w| — > oo, we easily 
see that the high-frequency limits of H and G are the 
same as for dielectric material, thus fijl 



lim — G(r,r',ui) = -8(r-r'). (A15) 
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To find the low-frequency limit of G, we note that the 
second term in Eq. I|A13|) is regular. To study the first 
term, we distinguish between two cases. 

(i) The first term of H in Eq. (|A1|) is transverse, 

— j"p x K(r,w)px = — p x n{r,u)f> x I" = 0, (A16) 

and therefore does not contribute to I^HlK It then 
follows that the same low-frequency behaviour as in the 
case of dielectric matter is observed, thus 0] 

UJ 2 

lim —G = -P[I^e(r,u) = 0)/"] -1 /", (A17) 

|u»|-»o c z 

i.e., because e(r, u = 0) ^ 0, 

lim — G(r, r', w) = M, My < oo. (A18) 

(ii) Equation 1A16I) is not valid, so that the first 
term of H in Eq. IjAlf) contributes to I" HI". Since 
(J,(t,u> — 0)^0 [thus k(t, uj = 0) being finite] , we find that 



lim G(r,r',w) = N, 



Nij < oo. 



(A19) 



APPENDIX B: COMMUTATION RELATIONS 

By using Eqs. (|22[1 . (|32|) . the commutation relations 
(|30|l and l|31|) . and the integral relation (|25fl. we derive 



[4(r^),4 T (r',c')] 
hw 2 



tteqc 



2 ImG«(r ) r' ) w)i(w-u') ) (Bl) 



[4(r,^),4(r',c')] - [^(r.wj^t^c')] = 0. (B2) 

From Eqs. Ip51)l. l]Blfl. and (|B2fl it is easily seen that the 
commutation relations l|35|l are valid. Moreover, we find 
that, on recalling that G*(r, r',w) = G(r, r', — u*), 



[e E i (r),M*')] 


= / d3 » 


"2*/i /■°° ( 


-/*. 


-7? / 

.TT '-oo 




°° dw W 2 
-oo W C 2 



dw -rImGijt(r, s,w) 



(B3) 



("P, principal part). Since the Green tensor is analytic in 
the upper half of the complex w-plane with the asymp- 
totic behaviour according to Eq. IjAlSjl . the frequency 
integral in Eq. I|B3|I can be evaluated by contour integra- 
tion along an infinitely small half-circle around w = 0, and 



along an infinitely large half-circle \uj\ — >oo. Taking into 
account that the Green tensor either has only longitudi- 
nal components in the limit \ui\ — ► 0, cf. Eq. (|A17|) . and 
hence uj/cGJ^-^O, or is well-behaved, cf. Eq. (|A19fl . we 
see that the integral along the infinitely small half-circle 
vanishes. Recalling Eq. (|A15|) . we then readily find 



[soE^Mr')} =iM±(r-r>). (B4) 
Since B(r) = V x A(r), from Eq. (|B4|) it follows that 

[e Ei(r), 4 (r')] = -ihe ljk d r k 5(v - r'), (B5) 

i.e., Eq. It is then not difficult to see that the 

commutation relation 1|B4|) implies 

&(r),Mr>)] = [<p(r),B i (r')]=0. (B6) 

To evaluate commutators involving the displacement 
field and the magnetic field, recall Eqs. I|16|) - i|19|) . Using 
the relations presented above, we derive 



and 



[A(r),A-(r')] =0, 
[H l (r),H J (r')]=0, 

[A(r),Mo4-( r 0] = -ihe ijk d r k 6(r - r'). 



(B7) 
(B8) 

(B9) 



Note that Eq. I|B9(I follows by using similar arguments as 
in the derivation of Eq. I|B4|) from Eq. (|B3|) . A similar 
calculation leads to 

[A(r),i j (r')]=iM i J(r-r'), (B10) 

[H^Ajir')] = [A(r),^(r')] = [^(r)^(r')] = 0. 

(Bll) 



By combining Eqs. {BUJ - ijETTjl with Eqs. l|B4|) and 
(|B6|) . it is not difficult to verify that polarization and 
magnetization commute with the introduced potentials 
as well as among themselves. 



APPENDIX C: HEISENBERG EQUATIONS OF 
MOTION 



By using the Hamiltonian (|41|) and recalling the def- 
initions of the medium-assisted field quantities in terms 
of the basic fields fxi(v,u>) the basic-field commutation 
relations (|30|l and l|31|l . the commutation relations that 
have been derived from them, and the standard commu- 
tation relations for the particle coordinates and canonical 
momenta, it is straightforward to prove that the theory 
yields both the correct Maxwell equations (|47|l and i|48|) 
and the correct Newtonian equation of motion 1)5 

Let us begin with the Maxwell equations. We derive 
on recalling Eqs. gUJ) and (O together with Eqs. 
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and 1)32)1 and the commutation relations 1)30)1 and (|31(l , The first term on the right-hand side of Eq. 1)C6)I is again 



B(r,t) = -[B(r,t),H] 



V x / duj — [A(r , u), H] + H.c. 

Jo l " 



■^fd 3 r Tdujhu; ]T [A(f a ),ft(r, W )f A (r,w)] 

"J J \ „, 



= - V x E(r) = -V x £(r), (CI) 

which is Eq. (|47|l . To derive the equation of motion for 
the displacement field, we have to consider several com- 
mutators according to 

5(r) = i[B(r,t),ff] 

= ^ Jd\' J^dujhu ^ [D(r),ft(r>)f A (r>)] 



;fcE^ [D(r),[p«-feA(r a )] : 



A— e,m 

iujq a k{r a ) = (?Q E i (f Q ). (C7) 



The second term gives rise to two terms, 



ig Q {f Q A(f Q )(g) V + V®A(f Q )f Q } (C8) 



a a 

The first commutator in Eq. I)C2)1 can easily be found 
by recalling the definitions of displacement and magnetic 
fields as 

^ / dV /°° dwftw E [D(r),fl(r',c)f A (r', W )] 



and 



-|E2^[ A ^)'[^-^ A ^) 



/ dwiwD(r,u)+H.c. = Vxi(r). (C3) = -§9a {fa V <g> A(f Q ) + A(f Q ) ® f Q v} , (C9) 

Jo 



Applying the commutation relation l)B10)l . and recalling 
the definition of the current density, we find that the 
second term on the right-hand side of Eq. ()C2)I can be 
written as 



and thus 



1 



1 



E PW' [P« - 9«A(f Q )] J = -ji(r), (C4) 

where the Newtonian equation of motion (|50|l has been 
used which follows directly from the Hamiltonian 1)41(1. 
Finally, standard commutation relations together with 
the definitions of the scalar potential, Eqs. 1)42)1 and 1)43 [I. 
and the current density, Eq. 1)491) together with Eq. (|50|l . 
lead to 



5«o 



r a x%)-%)xr Q . (CIO) 



£0 



By means of Eqs. H39|) and (|42|l one can see that the last 
two terms in Eq. (|C6|) can be rewritten as 



(C5) 

Inserting Eqs. I)C3)1 - l)C5)l into Eq. I)C2)1 . we arrive at 
Eq. 

In order to prove Eq. (|51H . we consider the equation 
= ^[p a - g a A(r a ), ill 



Jd 3 rj"dcjtiw ]T [A(r a ), (r, w)f A (r, w)] 



ifi E 2^ [P° _ <7aA(r Q ), [P/3 ~ ?/3 A (f^)] " 



i /" d 3 r [p Q)/ 5 A (r)^ A (r)] = -g Q V^ A (r Q ), (Cll) 
i 1 d 3 r [ Pai p A (r)^(r)] = <? q EU (r Q ). (C12) 



Inserting Eqs. iCT0|) - jHUt into Eq. |U6j and 

making use of Eq. 1)44( 1. we just arrive at Eq. 1)51)) . 



+ 



J d 3 r [p a ,p A (r)(£ A (r) 
i / d 3 r [p a ,p A (r)0(r)]. 



(C6) 
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